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ABSTRACT 

We present new X-ray and optical spectra of the old nova CP Pup (nova Pup 1942) ob- 
tained with Chandra and the CTIO 4m telescope. The X-ray spectrum reveals a multi- 
temperature optically thin plasma reaching a maximum temperature of 36Rg keV 
absorbed by local complex neutral material. The time resolved optical spectroscopy 
confirms the presence of the ~1.47 hr period, with cycle-to-cycle amplitude changes, as 
well as of an additional long term modulation which is suggestive either of a longer pe- 
riod or of non-Keplerian velocities in the emission line regions. These new observational 
facts add further support to CP Pup as a magnetic cataclysmic variable (mCV). We 
compare the mCV and the non-mCV scenarios and while we cannot conclude whether 
CP Pup is a long period system, all observational evidences point at an intermediate 
polar (IP) type CV. 
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1 INTRODUCTION 

CP Pup exploded in 1942 as one of the brightest and fastest 
(t 3 = 8 days, Schaefer & Collazzi 2010) of the Galactic clas- 
sical novae. The luminosity of the post-nova never returned 
to the pre-outburst quiescent value. The nova remained 
in a sort of standstill and, according to Schaefer & Col- 
lazzi (2010), CP Pup might still be more than five magni- 
tudes brighter than before the outburst. Thanks to its rel- 
atively bright “quiescent” magnitude and outburst charac- 
teristics, CP Pup has been widely studied as a rather “ex- 
treme” example of a nova, one that really helps constrain 
the parameters used in the classical nova (CN) models. The 
main findings were: 1 ) an unstable spectroscopic period of 
~1.47hr (Bianchini et al. 1985 and 2012, Duerbeck et al. 
1987, O’Donoghue et al. 1989, White et al. 1993); 2) a ~2 
to 11% longer (and unstable) photometric period (Warner 
1985, Diaz & Steiner 1991, Patterson & Warner 1998); 3) an 
X-ray spectrum similar to those of magnetic CV (mCV, Orio 
et al. 2009) , 4) a modulation of the X-ray flux at the spectro- 
scopic period (Orio et al. 2009, Balman et al. 1995); 5) dou- 
ble humped NIR light curve with period roughly matching 
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the photometric period (Szkody & Feinswog 1988). Taken 
together, these findings do not deliver a consistent picture of 
the old-nova (see details in Sect. 5), and, in particular, leave 
unresolved the debate upon the nature of the nova (whether 
magnetic or non-magnetic) and its mass (too small to be 
consistent with either classical nova and/or stellar evolution 
theories) . 

In this paper we present new Chandra X-ray data and 
time resolved optical spectroscopy which we collected with 
the aim of addressing those inconsistencies. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 The Chandra observations 

In the previous X-ray observation with XMM-Newton (Orio 
et al. 2009), the upper limit of the RGS grating (0.33- 
2.4 keV) left an uncertainty concerning the highest plasma 
temperature in the cooling flow model with which we fit- 
ted the hot plasma. This temperature indicates the depth 
of the gravitational potential well of the white dwarf (WD), 
and, therefore, a way to measure the WD mass. In order 
to better constrain the plasma temperature and the WD 
mass we observed CP Pup with Chandra ACIS-S camera 
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Table 1. The seven Candra exposures and their average count 
rate. Start and end times are in Terrestrial Time (TT). 


Obs.ID 

Date 

Start 

End 

Exposure 

(s) 

HEG 

counts/s 

MEG 

counts/s 

9967 

2009-09-30 

15:25 

21:29 

19814 

0.0196 

0.0356 

11990 

2009-10-02 

05:23 

19:20 

47263 

0.0207 

0.0360 

9966 

2009-10-11 

08:30 

16:01 

24729 

0.0198 

0.0367 

12005 

2009-10-15 

09:53 

19:41 

32504 

0.0237 

0.0430 

12012 

2009-10-26 

00:43 

06:53 

19715 

0.0246 

0.0482 

11991 

2009-10-28 

14:44 

20:52 

19714 

0.0252 

0.0449 

12013 

2009-10-31 

06:28 

09:45 

9187 

0.0225 

0.0445 


coupled with the High Energy Transmission Grating Spec- 
trometer (HETGS). The HETGS is a system of two sets 
of gratings, the HEG (high energy grating) and the MEG 
(medium energy grating). This allows high resolution spec- 
tra (with E/AE up to 1000) between 0.4 keV and 10.0 keV. 
CP Pup was observed between September 30 and October 
31 2009: the total observing time was 175,780 s, split in 7 
exposures because of scheduling constraints. Table 1 reports 
the log of observations together with the average count rate 
recorded during each exposure. 

2.2 The CTIO observations 

Spectroscopic observations were taken at the CTIO’s 4m 
telescope equipped with the, now retired, R-C spectrograph. 
The instrument setup was grating KPGL1-1 at a tilt angle 
of ~57° and slit width of ~1.2arcsec which, together, cov- 
ered the wavelength range 3500-6000 A at a resolution of 
~6 A x (i.e. ~320km/s). The exposures time was 60s and 
the duty cycle was ~ 1.8 min. A few times during each night 
the observing sequence was interrupted in order to take a 
wavelength calibration exposure. CP Pup was observed on 
Feb 6 and 7 2009 for about 5.9 and 4.7 consecutive hours, 
respectively. In total we collected 177 spectra during the 
first night and 125 during the second. The wavelength cali- 
bration exposures were used to calculate the dispersion so- 
lution of the spectra, however the wavelength “zero-point” 
was determined using the [O i] sky line at 5577 A, due to the 
spectrograph instability. Each spectrum was corrected also 
for the earth motion, in addition to the instrumental flexure. 

The data were reduced following the standard proce- 
dure and using the dedicated IRAF packages. All the spectra 
were flux calibrated using the spectrophotometric standard 
stars HR4468 & LTT4861 and LTT2415 for the first and the 
second night of observation, respectively. 


3 ANALYSIS OF THE CHANDRA DATA 

We reprocessed all the HETG data using chandrarepro in 
the CIAO v4.3 software package. Using the resulting pha2 
files, with standard source and background regions, we ex- 
tracted type I standard pha files by means of the ftool 
cmppha for the source and the background, and combined 
the ±l st orders into a single MEG file and a single HEG 
file, each with associated response and background (mak- 
ing use of the tool addascaspec). We then combined these 

1 This is the measured FWHM of the sky emission line. 


files from the 7 pointings into a single HEG spectrum and a 
single MEG spectrum. We focused on the ranges 0.4-6 keV 
(2.07-31.00 A) for the MEG and 0.7-8 keV (1.55-17.71 A) 
for the HEG. As the summed spectra have a modest sta- 
tistical quality, we produced two different versions of the 
spectra with different binning, using the ftool grppha. In 
particular, one was heavily re-binned so that each channel 
had at least 25 source photons; we fitted these files using 
X 2 statistic, which is well suited to evaluating the goodness 
of fit. We also used this version to produce the final plots 
in Fig.l. For the final fitting, done with XSPEC, vl2.6, we 
used a binning by a factor of 2, and the C statistic (Cash, 
1979), with the background-subtracted spectra. 

The spectra clearly contain K shell emission lines of var- 
ious species from Fe to O, as expected for X-ray emission 
from a multi-temperature plasma (Fig. 1). The strongest 
lines are those of He-like Fe at 1.86 A, He-like Ar at 3.96 A, 
H and He-like S at 4.73 A, and 5.06 A, respectively, H and 
He-like Si at 6.17 A, and 6.67 A, respectively, H like Mg at 
8.43 A, and H and He like Ne at 12.16 A, and 13.47 A, re- 
spectively. Hence, following Mukai et al. (2003), we modeled 
the spectra using a cooling flow model, and fixed the mini- 
mum temperature to the lowest value available in the model, 
kT=0.0808 keV. This is effectively zero, as cooler plasma 
does not contribute with significant emission in the HETG 
band. 

In fitting the models with XSPEC, we started by using 
a phabs (mkcf low+ga) model, with the Gaussian component 
representing the fluorescent Fe line at 6.4keV (1.86A). How- 
ever, the fit was poor so that we applied a series of modi- 
fications to the model. We first added a partial covering 
absorber improving the fit significantly. Second, since the 
line strengths deviated from the predictions of a solar abun- 
dance cooling flow model, we used the variable abundance 
version of the cooling flow model (vmcf low) and allowed the 
abundances of O, Ne, Mg, Si, S, and Fe to vary (we tied the 
Ni abundances to Fe). Finally, we applied the gsmooth con- 
volution model to the vmcflow component with its second 
parameter set to 1.0 (constant velocity width) to broaden 
the model line widths to match the observations. 

Applying the best-fit model from the C statistic fit 
to the heavily grouped spectra, we obtained a reduced \ 2 
of 1.1048 for 416 degrees of freedom (and 1.0870 for 413 
degrees of freedom above 0.5 keV). Many of the fit pa- 
rameters (27 in total 2 ) are poorly constrained. In partic- 
ular, we cannot derive detailed elemental abundance but 
rather get the abundance pattern, which is suggestive of 
Ne abundance about twice solar (i.e. 2.2F0.6). We also de- 
rive the column density for the simple absorber and the 
partial covering absorbers: Nff=0.16(±0.03) x 10 22 cm -2 and 
3.0(±0.66) x 10 22 cm -2 , respectively; with a covering fraction 
of 45% for the partial covering absorber. We notice that 
the H column density measured by LAB and GLASS from 
the 21cm transition (Kalberla et al. 2005, Kalberla et al. 
2010) in the direction of CP Pup is 0. 672-0. 673x 10 22 cm' 2 , 
i.e. significantly smaller than our total column density es- 
timated by modeling of the Chandra spectra. This is an 


2 They include the elemental abundances for 14 atoms, minimum 
and maximum temperature, column density of the simple and 
partial absorber, redshift, normalization/m, etc. 
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Figure 1. Top panel: MEG data and their best fit model (red line) & HEG data and their best fit model (blue line). Bottom panel: 
data/model ratio, with the same color convention as above. Note that the wavelength scale is logarithmic. 


independent demonstration that we are witnessing local ab- 
sorption in CP Pup, consistently with a complex local ab- 
sorber. The model (absorbed) flux is 4. Ox 10~ 12 erg/cm 2 /s 
in the range 2-10 keV and 4.8x 10 -12 erg/cm 2 /s in the range 
0.4-10keV. This corresponds to 4. 5x 10 -12 erg/cm 2 /s and 
6. 9x 10 -12 erg/cm 2 /s, respectively, once corrected for ab- 
sorption, and to a luminosity of 1.4e33 erg/s (2-10keV) and 
2.1e33erg/s (0.4-10 keV) for an assumed distance of 1.6 kpc 
(Williams 1982). For the same distance, the normalization 
constant in the cooling flow model indicates an accretion 
rate of 4.3x 10~ 10 Mg/yr (best fit) and a 90% confidence 
range 3.3-7.3X 10~ 10 Mg/yr. The strong lines appear broad- 
ened at the 200 km/s level. Finally, we also derive our main 
parameter of interest, the maximum temperature, which 
we use to constrain the WD mass and for which we ob- 
tain kT max =36.5 keV, with a 90% confidence range of 20.2- 
55.7 keV. 

The presence of the complex absorber is the hallmark of 
magnetic CV’s (Ramsay et al. 2008). In addition, the excess 
of counts below 0.5 keV could indicate a soft component aris- 
ing from the heated white dwarf (WD) surface. Hence, we 


computed the WD mass assuming that CP Pup is indeed a 
magnetic CV and obtained Mi=0.801q 23 Mq. Note that for 
maximum temperatures as high as ^36keV, small changes 
in kT max only lead to subtle difference in the continuum 
in the HETG range, so we must assume that we modeled 
the complex absorber sufficiently accurately. However, the 
best handle in the HETG data for kT max is the emission 
line strength. Plasmas with kT higher than ~20keV adds 
continuum flux but little line emission in the HETG band. 
Thus our constraint comes from the observed line equiva- 
lent width, and the large error range is partly due to the 
fact that we allowed the abundances of several elements to 
vary. 

Finally we have analyzed the X-ray data searching for 
possible extended emissions originating from the nova shell 
similarly to, for example, GKPer (Balman 2005). We were 
unable to detect any X-ray extended emission and our sur- 
face brightness profile computed on the spectra perfectly 
matches the simulated PSF, once the constant background 
level has been taken into account. 
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Figure 2. Left panels: a representative spectrum of our time series. It has been obtained by averaging two spectra taken toward the end 
of the first night. Note that the spikes at 4130 4517 5018 5065 5395 and 5665 A are all bad pixels. The inset compare the H (3 line profile 
(solid line) with that of the sky emission lines OiA5577 (dashed line). Right panel: example of profile variation shown by the Chi and 
Hen emission lines. The four spectra are separated by about 25-50 minutes each other. 



8.6 8.65 8.7 8.75 8.8 9.55 9.6 9.65 9.7 9.75 

(HJD^2454860.) 


Figure 3. Top panels: radial velocity curves of the emission lines H/3 (black solid circles) and H7 (blue circles) during the Feb 6 (left) and 
Feb 7 (right) observing runs. Bottom panels: radial velocity curve of the emission lines HeiiA4686 (blue solid circles) and the Cin+N 111 
blend (red circles) during the same nights. The y range is the same in all the panels. 
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4 DATA ANALYSIS: THE OPTICAL SPECTRA 

We measured the wavelength of the strongest emission lines 
in each spectrum, i.e. the Balmer lines H/3 and H7 as well 
as the Hen A4686 and C IIl(l). Given the low resolution and 
the signal-to-noise ratio (SNR) of the spectra we could not 
distinguish a broad and a narrow component of the emission 
lines (e.g. Bianchini et al. 2012, O’Donoghue et al. 1989) 
and fit each emission line with a single Gaussian profile, 
tracking just the dominant component. We note however, 
that the line profile, particularly in the case of the C ill and 
Hell emissions, occasionally deviates from Gaussian profile 
both because of intrinsic variations and poor SNR. Fig.2 (left 
panel) shows the average of two spectra taken at the end of 
the first night; while the figure’s inset, compares the profile 
of the H/3 emission (solid line) with that of the sky line 
OlA5577 (dashed line). The right panel of Fig.2, instead, 
shows the variability displayed by the Cm and the Hell 
emission lines. 

The radial velocity (RV) curve of each line is plot in 
Fig. 3. It shows that: 

(i) there is an overall drift during the first night and an 
offset of the average radial velocity between the two nights. 

(ii) the amplitude of the ~ 1.47 hr period may vary from 
cycle to cycle. 

(iii) there is a large scatter of the data points with occa- 
sionally deviating points which only in part can be explained 
by poor SNR and/or fitting. 

Point (i) and, in particular, the trend of the first night 
suggest a variability on timescale longer than our observing 
run. A Fourier analysis combining the whole data set of two 
nights gives a power spectrum that peaks at the well known 
period at ~16.25 cycles/day (~1.47hr), as well as a substan- 
tial power at low frequencies around ~2.38 cycles/day (cor- 
responding to ~10.1 hr). A formal fit to the H/3 radial veloc- 
ity curve gives 9.8 hr (with the uncertainty of 1 d, including 
all the aliases), which has to be regarded as a lower limit 
to a possible periodic variation. Although there is certain 
evidence of long term variability, we are unable to establish 
whether it is periodic or not due to the time gap between 
the two nights and the time coverage of our observing run. 

Points (ii) and (iii) suggest that the line forming region 
might be varying in location, geometry and/or optical depth. 

Phasing the spectra on the ephemeris: 
T 0 =2452868. 582093+nx0. 0614466 HJD, which we de- 
rived from the best fit for the H/3 RV curve, we observe 
that (see Fig. 4 top panels) different ions have different 
K-amplitudes, different average-RV and same red-to-blue 
crossing phase, suggesting a structured line forming region. 
The line dispersion velocity (FWHM) is fairly large in 
the range 1000-2000 km/s and modulated across the high 
frequency period (see Fig. 4, bottom panels). The modu- 
lation varies in amplitude and shape with the ions, again 
consistently with a structured line emitting region and 
possible optical depth effects. 

We remark that we did not observe modulations of the 
lines equivalent width (EW), nor correlation between, e.g., 
the Hell and the H/3 EWs. This possibly indicates that the 
line strength is larger when a more collimated area is into 
view and/or the presence of more than one line forming 
region. We also note that the C ill velocities and variable 
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Figure 4. Top panels: radial velocity curve of the emission lines 
H/3 and H 7 (left) and He II and C III (right - same color code as in 
Fig. 3). Bottom panels: modulation of the emission lines FWHM 
with the ~1.47hr period. Note that while H/3 and H 7 FWHM 
present an almost sinusoidal modulation, the He II modulation is 
more similar to a sawtooth profile. The C III FWHM is much less 
modulated and “noisy”, possibly indicating a variable contribu- 
tion/blend with the N III emission line. 

profile are suggestive of a significant contribution from the 
N IIIA4640 line of ~25%. The lack of significant modulation 
of its FWHM might indicate that the N ill fractional contri- 
bution is phase dependent. 

Finally we plot in Fig. 5 the trailed spectra of the var- 
ious emission lines folded on the 1.47 hr period. In addition 
to the “standard” trailed spectra we produced “difference 
trailed spectra” by subtracting CP Pup average spectrum 
to each spectral bin. Difference trailed spectra enhance the 
line asymmetries and their motion. The trailed spectra re- 
veal a “helix” pattern of the emission lines across the 1.47 hr 
period which became suggestive, again, of a structured or a 
multicomponent line forming region particularly in the case 
of the He 11 line. 


5 DISCUSSION 

In Section 1 we mentioned that the observations collected at 
all wavelength (X to NIR) for the quiescent CP Pup provide 
inconsistent results. On one side the X-ray observations al- 
ways favored a magnetic nature of the nova: the X-ray flux 
is modulated at the high frequency period (Balman et al. 
1995, Orio et al. 2009) suggesting that it corresponds to the 
WD spin period. In addition, the observed X-ray spectrum 
is best modeled by high temperature optically thin plasma 
that is commonly found in magnetic systems. On the other 
hand, polarimetry observations (Cropper 1986) did not re- 
veal polarized light in CP Pup disfavoring the polar inter- 
pretation, though not the IP scenario. However, Cropper’s 
(1986) observations cannot be conclusive due to their short 
time coverage (18 min). Time resolved optical observations 
(Bianchini et al. 1985, 2012, Warner 1985, Duerbeck et al. 
1987, O’Donoghue et al. 1989, Barrera & Vogt 1989, Diaz & 
Steiner 1991, White et al. 1993, Patterson & Warner 1998) 
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Figure 5. Trailed spectra (left panels) and “difference trailed spectra” (right panels) of the emission lines H/3, Hq, Hell + Cm and He I, 
from top to bottom. Phase 0 corresponds to the time for the red-to-blue crossing in the high-frequency modulation. 


detected unstable periods <2 hr. In particular, spectroscopic 
periods in the range 0.06115-0.06148 days, were interpreted 
as the binary orbital period, while photometric periods in 
the range 0.06138-0.0614 days were interpreted as super- 
humps, as commonly observed in short-orbital period sys- 
tems accreting at high rates (Patterson & Warner 1998). In 
view of these periods CP Pup was recognized as the first CN 
below the period gap. However, the WD mass value deter- 


mined from optical RV curves is uncomfortably low, in the 
range 0.02^Mi <0.6 Mq depending on the author and the 
adopted orbital inclination (White et al. 1993, O’Donoghue 
et al. 1989, Duerbeck et al. 1987, Barrera & Vogt 1989, Bian- 
chini et al. 2012). Only Bianchini et al. (2012) were able to 
force an analytical solution with Mi ~1 Mq by including a 
peculiar hot spot contribution to the line flux. However, no 
past, nor present observations are consistent with the con- 
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figuration of accretion disk + hot spot. WD masses <0.6 Mq 
are hardly compatible with the classical nova (CN) theory 
and, in particular, with the fast super-Eddington outbursts 
of the type displayed by CP Pup. These outbursts require 
more massive WDs, typically^O. 8-1.0 Mq (e.g. Starrfield et 
al. 2012a, 2012b, and reference therein). Finally the NIR 
J-band light curve shows a double humped modulation at 
about the photometric/spectroscopic period, which, if inter- 
preted as ellipsoidal variations, implies a secondary star that 
is not compatible with a short orbital period system below 
the period gap and at the distance 3 of 1.6 kpc of the nova. 
CP Pup NIR. magnitudes are also incompatible with the bi- 
nary evolutionary phase at the measured distance, being far 
too bright. 

Our new Chandra and CTIO observations support the 
magnetic nature of CP Pup. The Chandra data show ev- 
idence of a multi-temperature plasma which is absorbed 
by dense material partially covering the X-ray source and 
temperatures up to kT=36.5keV. High plasma tempera- 
tures and complex local absorbing medium are characteris- 
tic of mCVs (Ramsay et al. 2008). Non-mCVs do not suffer 
strong local absorptions except for high orbital inclination 
systems (e.g. the eclipsing system V893Sco, Mukai et al. 
2009). Note that only the Chandra HETG has the sensi- 
tivity, band-pass and spectral resolution that allow to as- 
sess the presence of a complex absorber. Note also that the 
WD mass can be derived from the maximum plasma tem- 
perature, kT m ax = 36.5 keV, either regarding it as shock 
temperature onto a magnetized WD (i.e. free-fall), or as 
a temperature produced by accretion from a disk onto a 
non- magnetized WD (i.e. Keplerian velocities). In the lat- 
ter case the kT max temperature would be reached by a 
more massive WD: Mi = l.llg'ggMg (as opposite to the 
Mi = 0.8 ^q'23 Mq in the magnetic WD hypothesis). Both 
values, together with their uncertainties fall within the typ- 
ical CV primary masses (Zorotovic et al. 2011) and point to 
a higher WD mass than previously estimated which is more 
consistent with the CN theory. 

Our sequence of optical spectra imposes some strong 
constraint on the location and geometry of the line form- 
ing region suggesting a composite line forming region with a 
temperature gradient off the orbital plane. In addition, the 
observed 7 -drift implies either that the RV are affected by 
non orbital mechanisms, or the existence of a longer period 
(yet to be confirmed). Note that the 7 -drift is not just a one 
night occurrence as it appear also in the re-examination of 
the older data sets (Bianchini et al. 1985, 2012). It is intrigu- 
ing to notice, however, that we mostly detect positive drifts. 
Whether this is due to our specific sampling or is intrin- 
sic to the line emitting source cannot be concluded at this 
stage, nor we can say whether the phenomenon is strictly 
periodical. Anyway, the lines must form in the binary ac- 
cretion component. They cannot, for example, form on the 
irradiated face of the secondary star for their FWHM are 
too large to be consistent with such a hypothesis and the 
emitting secondary cannot explain the 7 -drift. In addition, 
the broad and narrow components identified by O’Donogue 
et al. (1989) and Bianchini et al. (2012) and which we do 
not resolve, appear to be in phase arguing, too, against sec- 


3 Calculated from the nebular parallax (Williams 1982). 


ondary star irradiation. The accretion component can either 
be an accretion disk in a non-mCV or the magnetically fun- 
neled gas in an mCV. Below we discuss the two scenarios 
separately in the attempt to identify which one apply to 
CP Pup. 

non-mCV 

The classic disk + hot-spot configuration explains the differ- 
ent 7 velocities from the different lines, but not their match- 
ing phases, nor the 7 -drift. The hot-spot typically shows 
temperature gradients along the stream trajectory on the 
orbital plane imposing significant phase offsets among dif- 
ferent lines (e.g. Mason et al. 2000 and reference therein), 
contrarily to CP Pup observations. It is possible to observe 
different 7 velocities for a same line at different epochs. This 
is explained with long term variation of the geometry and/or 
velocity field of the line forming region, possibly associated 
to changes in the mass transfer rate. However, “real time 
variations” as those observed in CP Pup, if due to sudden 
mass transfer rate changes, should be associated to signifi- 
cant line intensity and profile variations as observed during 
flares or episodes of intermittent mass transfer (e.g. Mason 
et al. 2007, Rodriguez-Gil et al. 2012). We observe only small 
scale variations which are consistent with the period modu- 
lation and the night conditions. An eccentric and precessing 
accretion disk could explain the 7 -drift and, in fact, produce 
a sinusoidal modulation of the 7 velocity. The precessing 
disk would be consistent with Patterson & Warner (1998) 
super-hump proposal to explain the light curve modulations, 
though their instabilities (the “anomalous super-hump”, as 
defined by Patterson & Warner 1998), remain quite unex- 
plained. We note that this model is no longer valid should 
future observations prove that the 7 modulation follow a 
“sawtooth” pattern instead of a sinusoidal one. The exact 
location of the line forming region on the disk + hot-spot, 
however, remain unclear as there is no obvious configuration 
for having smaller K-amplitude in case of higher ionization 
potential ions. The observed FWHM are quite large for a 
low orbital inclination system and are hard to explain as 
are their modulation. Our trailed spectra do not match the 
pattern displayed by disk + hot-spot or hot-spot emission. 

mCV 

First of all, the presence of strong Hell emission lines is 
a strong indicator of the magnetic nature of the CV (van 
Paradijs & Verbunt 1984). Second, in a mCV the bulk of the 
line emitting region is located in the magnetically confined 
accretion flow toward the WD magnetic poles (Schwope 
et al. 2000). This configuration naturally imposes an off- 
orbital-plane/vertical structure characterized by tempera- 
ture and velocity gradients and non-Keplerian velocities. 
The various emission lines are displaced along the curtain 
and therefore their RV curves have matching phases but 
different 7 velocities. The higher ionization potential en- 
ergy transitions form closer to the magnetic pole and the 
K-amplitude distribution results from a projection effect. 
Similarly the line FWHMs and their modulation are im- 
posed by the velocity field along the magnetic field lines 
and its projection on the line of sight and consistent with 
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those expected in IP (Warner 1995 and reference therein, 
Ferrario et al. 1993). The helix pattern shown by the trailed 
spectra suggest multiple components of the type observed 
in mCV and IP in particular (see e.g. Hcllier 1999). Finally, 
the 7 velocity drift suggests two different possible scenarios: 

1) the modulation is sinusoidal and induced by the orbital 
motion of the WD around the binary center of mass. I.e. the 
~ 1.47 hr period is the WD spin and the true orbital period 
is several hours long (similarly to RXJ154814. 5-452845, see, 
e.g. de Martino et al. 2006) and has yet to be determined; 
while the photometric period is a beat. 

2) the 7 modulation is asymmetric and sawtootli-like and 
the line emitting region is affected by cyclical disturbances. 
We note at this point that neither the WD spin, nor the 
binary orbital period are expected to be unstable. However, 
they can appear as such if their probe is unstable. The accre- 
tion curtain (or threading region) is fixed in the binary frame 
if the system is synchronous; it is not otherwise. The curtain 
will have constant geometry and emissivity if it is fed by a 
symmetric accretion disk/ring, but will vary in case of asym- 
metries within the disk (either because the disk is missing or 
because the accretion stream is deeply penetrating into it). 
Hence the observer will view a line emitting region which 
is modulated with the WD spin and its beat. We note to 
this regard that different accretion geometries are possible 
in IPs, depending on the magnetic field strength and the 
ratio between the spin and the orbital period (Warner 1995 
and reference therein, Norton et al. 2004). We also note that 
the magnetic system BY Cam shares a number of similarities 
with CP Pup. It has a record of multiple and non repeatable 
periods with no optical nor polarimetry observations match- 
ing either the spin or the orbital period and depending on 
the sampling (e.g. Silber et al. 1997 and reference therein, 
Mouchet et al. 1997). Both light and RV curves show a phase 
shift and an abrupt reset to 0 which has been explained with 
the change in the magnetic geometry seen by the accretion 
stream as the WD rotates relative to the secondary and the 
switching of accreting magnetic pole 4 . It is not difficult to 
imagine that similar configuration could cause a drift in 7 
velocity beside the phase. 

Independently on the value of the orbital period, the 
double humped NIR. modulation might be explained with 
irradiation (of the secondary star or the accretion compo- 
nent, e.g. Penning 1985) by two accreting magnetic poles. 
Note that CP Pup is a low orbital inclination binary, hence, 
we might very well miss the direct view of the lower magnetic 
pole (or see its periodic eclipse while the upper pole is always 
visible), but still see its effects on the secondary (/accretion 
component). Strictly speaking the double humped modula- 
tion in this case should result in uneven maxima, however 
this is hard to say from Szkody & Feinswog (1988) plot due 
to the paucity of data points and the lack of error bars. 
Most important those data cannot be conclusive on irradia- 
tion (nor on ellipsoidal variations), for their low amplitude 
modulation might be not significant and their phasing with 

4 Mouchet et al. (1997) invoke a precessing WD, but in facts 
this is unnecessary. At the same time a similar effect could be 
obtained without invoking the change of accreting pole but the 
reconnection of the magnetic field line from the opposite side once 
the WD magnetic pole has swept over 180° with respect to the 
LI point. 


respect to optical photometry is unknown. Note also that 
lack of detection of secondary star features does not argue 
against CP Pup as a possible long orbital period system. On 
one hand, the system spectroscopic observations have always 
been limited to the Hell H/3 wavelength range, which is 
devoid of secondary star absorptions, and to low spectral 
resolutions, which smear out any possible weak absorption 
feature. Our CTIO spectra do not have the sufficient spec- 
tral resolution and SNR, too. On the other hand, the chance 
for detecting the secondary depends also on the tempera- 
ture and size of the various emission sources into play and 
it might well be that the whole accretion component com- 
pletely outshine a relatively cold secondary star. 

In summary, while we cannot be conclusive on the 
strictly-periodic vs semi-periodic nature of the observed 7- 
drift and its modulation, comparing the existing observa- 
tional characteristics of CP Pup within the mCV and non- 
mCV scenarios, the nova IP nature best fit them all by far. 
This is shown in Table 2 where we list CP Pup observational 
facts and how each scenario succeed or fail in explaining 
those consistently. 

For completeness we should also consider that the pos- 
sible existence of a longer period could imply that CP Pup is 
in fact a triple system, i.e. a CV below the period gap “orbit- 
ing around” a third body. In this case the hypothesis of an 
eccentric orbit around the third body would be needed if the 
asymmetry of the long term RV modulation with slow/long- 
lasting velocity increase and short /rapid velocity decreases, 
is confirmed. This scenario however is somewhat cumber- 
some and ad hoc and would open questions about its prob- 
ability and stability. It would probably requires formation 
within a globular cluster with subsequent ejection, an evo- 
lutionary scenario which, beside being very unlikely, needs 
detailed simulation (see, e.g. Anosova 1996, Brasser 2002) 
that are not possible with the current set of information in 
hand. Hence we do not pursue this scenario any longer and 
note that 1) the only CV that has been suggested to be 
part of a triple system is FSAur with a period around the 
third body of few hundreds of days (Chavez et al. 2012); 
2) placing a non-mCV in a triple system we will face the 
same contradictions and problem we are discussing for the 
isolated non-mCV. 


6 SUMMARY AND CONCLUSIONS 

In this paper we presented and analyzed Chandra/HETG 
X-ray spectra and high cadence, time resolved optical spec- 
tra obtained with the Blanco telescope at CTIO. The two 
data sets together produce interesting results, demonstrat- 
ing that CP Pup is very likely to be an IP system. The X-ray 
spectral characteristics are consistent with those of known 
IP’s. Assuming that CP Pup hosts a magnetic WD we fitted 
the X-ray spectrum with a cooling-flow model and a com- 
plex absorber, and from the best fit we derived a mass of 
0.8 Mq for the WD, which is in reasonable agreement with 
the typical WD masses in CV and with the prediction from 
the CN outburst theory. The IP hypothesis thus solves the 
problem of the non-realistic low mass derived assuming that 
accretion occurs only through a disk. 

Although in Orio et al. (2009) some of us derived an 
upper limit of 1.6 x 10 -10 Mq year -1 for the mass accretion 
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Table 2. Comparison of the possible scenarios applicable to CP Pup: the non magnetic CV, the IP (or asynchronous polar), both as 
long or short orbital period system. The first column lists the observed characteristics and for each model we report whether it fits the 
observation (Y) or not (N) or whether it fits them just partially (~). We also describe how the model explains/fits the observation or 
why it does not, as appropriate. 


OBSERVATION 

non-mCV 
short period 

mCV/IP 

long orbital period 

mCV/IP 

short orbital period 

strong HeiiA4686 

[N]: quiescent (low m) non-mCV 
typically do not show Hell emission 

[Y] : strong He II emission line is a 
characteristic of mCVs 

[Y]: strong Hell emission line is a 
characteristic of mCVs 

long term trend 

[Y] : the drift in 7 is induced by a 
precessing eccentric disk. 

[Y] : the observed changes in the 
7 velocity corresponds to the true 
(still undetermined) orbital period 
of CP Pup; the ~1.4 hr period is the 
WD spin. 

[Y]: the 7 drift is induced by the 
varying geometry of the magnetic 
field lines as seen by the accretion 
stream due to the WD rotation. 

NIR light curve 

[N]: ellipsoidal variations due to 
the elongated secondary star. El- 
lipsoidal variations are inconsistent 
with both the distance and the evo- 
lutionary phase of the binary. 

[Y] : the spinning WD magnetic 
poles outshine twice per WD ro- 
tation the secondary star (/accre- 
tion component) which reflects and 
reprocesses the X-ray similarly to 
DQ Her. 

[Y]: as for the long orbital period 
mCV scenario. 

Difference between spectro- 
scopic and photometric period 
(by a few %) 

[Y]: the photometric period is the 
super- hump period, while the spec- 
troscopic period is the orbital one. 

[Y]: beat between spin and orbital 
period. 

[Y]: phase jitter and unstable peri- 
ods as observed in BY Cam. 

dynamical WD mass 

[N] : the primary masses derived 
from dynamical solutions assuming 
that the optical emission lines trace 
the binary orbital period are always 
too low to be consistent with stellar 
evolutionary theory and/or the CN 
outburst theory. 

[Y]: reasonable masses could be 
found with the longer orbital pe- 
riod. 

[Y]: the emission line region simply 
does not trace the orbital motion 
of the binary and cannot be used 
for the computation of dynamical 
masses. 

X-ray spectrum 

[N]: the observed spectral energy 
distribution and high temperature 
match typical values of mCVs. 

[Y]: the spectral characteristic are 
those of an mCV 

[Y]: as for the long orbital period 
mCV scenario. 

presence of complex absorber 
in the high energy X-ray 

[N] : this is not typical of non-mCV 
and, in particular, low orbital incli- 
nation non-mCVs 

[Y]: this is an hallmark of mCVs 

[Y]: this is an hallmark of mCVs 

X-ray light curve 

[N] : in non-mCV disk systems the 
X-ray emission is from the bound- 
ary layer which is supposedly sym- 
metric producing no modulations. 

[Y]: the X-ray modulation is con- 
sistent with mCVs. The modulation 
result from the occultation of one 
magnetic pole. 

[Y]: as for the long orbital period 
scenario. 

modulation of the high fre- 
quency period (both photomet- 
ric and spectroscopic) 

[~]: anomalous super- humps 

[Y]: phase jitter and period insta- 
bility induced by asynchronism be- 
tween the WD spin and the orbital 
period, similarly to BY Cam. 

[Y]: as for the long orbital period 
scenario. 

modulated amplitude (spectro- 
scopic and photometric) 

[~] : anomalous super-humps 

[Y] : variable magnetic field geom- 
etry as viewed from the secondary 
star and the accretion disk, due to 
the WD rotation and chaotic accre- 
tion. 

[Y]: as for the long orbital period 
scenario. 

undetected circular polariza- 
tion 

[Y] : non-mCV do not show polar- 
ized light 

[Y]: IP polarization is typically 

masked by the disk/accretion com- 
ponent emission 1 

[Y]: as for the long orbital period 
scenario. 

undetected secondary star fea- 
tures 

[Y] : the secondary is a cool late 
type main sequence star as ex- 
pected for short orbital period 
systems. Note that this is “self- 
contradictory” with the explanation 
provided for the non-mCV scenario 
on the NIR light curve point. 

[~]: the secondary star should be 
relatively large and possibly com- 
peting in size (but not in tempera- 
ture) with the accretion component. 
It should be searched for at NIR 
wavelengths. 

[Y]: as for the non-mCV scenario. 


1 Note that the non-detection of polarized light favor the IP scenario more than the asynchronous polar one. 

^ Note that a long orbital period does not automatically implies that the secondary is visible in the R, I or NIR bands: high mass transfer rate and 

large/luminous accretion disk can mask the secondary star absorptions. 


rate m onto the CP Pup WD, the Chandra HETG spectrum 
indicates an m upper limit of 3.3 x ICC 10 Mq year -1 within 
a 90% confidence level, which may be too high to develop a 
very luminous TNR (Yaron et al. 2005) and leaves open the 
possibility of a variable rn before and after the outburst. 

Like the X-ray data, also the sequence of optical spec- 


tra, the emission line characteristics and their radial veloc- 
ities are much better explained with the IP scenario rather 
than with a non-mCV. We also observed a drift of the sys- 
tem 7 velocity, which is either suggestive the presence of a 
cyclical non Keplerian mechanism or of a longer period. We 
speculate that this “hidden” period may be the true orbital 
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one, rather than the spectroscopic period that may be as- 
sociated with the WD rotation. While with the current set 
of data in hand we cannot be conclusive about the origin 
of the observed 7 velocity drift, we conclude that we found 
very strong evidence that CP Pup is an IP. 

New optical observations are needed in order to solve 
the pending issues. High-cadence time resolved spectroscopy 
for several consecutive nights would allow to determine the 
period and the morphology of the 7 velocity modulation. 
The number of nights we requested at the Blanco was 3, 
but we lost one due to weather, and it turns out that a run 
consisting of more nights (>4) is necessary. A simultaneous 
runs of broad band, multicolor photometry done with an- 
other telescope should allow to assess correlations between 
light and RV curves, to observe whether the continuum and 
the line emission originate in the same region, and determine 
how many different sources of optical exist. High resolution 
spectra are desirable to properly identify and trace the var- 
ious emitting components. A long monitoring photometric 
program should be undertaken to best identify the beat peri- 
ods into play. While we are aware that for this ideal dataset 
many nights of observations with telescopes of different sizes 
and characteristics would be necessary, we do think that CP 
Pup is such an emblematic object for the nova theory that 
obtaining these results would be very rewarding for a full 
understanding of novae. 
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